Spatially-Resolved Temperature Measurements of Electrically-Heated Carbon Nanotubes 
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Spatially-resolved Raman spectra of individual pristine suspended carbon nanotubes are observed under elec- 
trical heating. The Raman G+ and G_ bands show unequal temperature profiles. The preferential heating is 
more pronounced in short nanotubes (2 /^m) than in long nanotubes (5 ^m). These results are understood in 
terms of the decay and thermalization of non-equilibrium phonons, revealing the mechanism of thermal trans- 
port in these devices. The measurements also enable a direct estimate of thermal contact resistances and the 
spatial variation of thermal conductivity. 



The temperature of a macroscopic solid is manifest in the 
energy of its lattice vibrations, or phonons. In nanostructures, 
this definition of temperature can break down if the phonon 
population is driven out of equilibrium by an electrical cur- 
rent. When that is the case, some phonon modes can have 
higher effective temperatures than the rest of the lattice de- 
pending on their coupling with the current. In the past, several 
interesting observations have been made in electrically-heated 
carbon nanotubes such as current saturation [l], F- and K- 
point optical phonon scattering ||2l O, ballistic phonon prop- 
agation I?] and electrical breakdown ||5] . In freely suspended 
carbon nanotubes, where phonon relaxation is hindered due 
to the absence of a substrate acting as a thermal sink, striking 
negative differential conductance (NDC) at high electric fields 
has been observed |6 1 and explained using non-equilibrium or 
"hot" phonons |6 7 |. 

These high-field properties are particularly relevant for ap- 
plications of carbon nanotubes as field-effect transistors and 
interconnects towards the miniaturization of electronics. Ad- 
ditional insight into these properties can be gained through 
simultaneous optical and transport studies. In particular, be- 
cause nanotubes are one-dimensional structures with a huge 
aspect-ratio, these phenomena can vary spatially, necessitat- 
ing a local probe of temperature in order to fully understand 
thermal transport in nanotubes. Previously, scanning force mi- 
croscopy |8 1 and local melting of nano-particles |9 1 have been 
used to extract local temperatures of multi-walled nanotubes 
under high bias. However, these contact-based techniques in- 
volve temperature drops at the measurement interface, ham- 
pering their ability to accurately measure temperature. 

Raman spectroscopy is a powerful, non-contact method of 
probing phonons in nanotubes |10|. This technique enables 
one to probe the F-point longitudinal optical (LO), transverse 
optical (TO) and the radial breathing mode (RBM) phonons, 
among others. Recently, hot phonons in nanotubes have been 
directly observed using Raman spectroscopy in conjunction 
with electrical transport 111 11 1121 . using the temperature re- 
sponse of the Raman G band in nanotubes 113, ,14, ,151. Thus, 
spatial investigation of electrical heating using Raman spec- 
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FIG. 1: Raman spectra of the center of device Dl at two bias volt- 
ages. Left inset: /-V characteristics of Dl. Right inset: SEM view- 
graph of a typical device. Scale bar corresponds to 2 /im. 



troscopy forms the motivation for our experiment. 

Individual single-walled carbon nanotube devices were 
grown on top of Pt leads, as reported in detail previously 1 16 |. 
Devices used in this work were grown using ethanol as the 
carbon feedstock, which has been shown to yield long and 
low-disorder nanotubes |17|. Figure 1 (right inset) shows a 
scanning electron microscope image of a 5 /im long nanotube 
device. Raman spectra were measured in a Renishaw In Via 
spectrometer with Spectra-Physics 532 run solid state and 785 
nm Ti sapphire lasers. Typical integration times were 60-120 
seconds. An Ithaco current preamplifier was used to measure 
the current through the nanotube. All measurements were per- 
formed at room temperature, and in an argon environment to 
prevent burn-out of devices at high bias voltages. Figure 1 
(left inset) shows the current-voltage (/-V) characteristics of a 
typical quasi-metallic device Dl (5 ^va long). Note the NDC, 
characteristic of suspended devices, above Vsd ~ 0.6 V in Dl. 

Figure 1 shows Raman spectra taken at the center of Dl us- 
ing a 532nm laser at two bias voltages, 14^=0 V and 1.2 V. 
Note the narrow G+ and broad G_ peaks that are characteris- 
tic of the TO and LO phonon modes, respectively, of metallic 
nanotubes. The two peaks downshift unequally in energy on 
application of a 1.2 V bias voltage. A low enough laser power 
was used such that the laser itself did not cause any downshift 
or heating. Note that the defect-induced D band peak in the 
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FIG. 2: (Color online) (a) Spatial profile of Raman shifts of G+ and 
G_ bands for device Dl at V and 1.2 V bias, (b) Temperature 
profiles of G+ and G_ bands. Inset: G+ and G- shifts for Dl in a 
temperature-controlled stage. 



Raman spectrum, occurring in typical nanotubes around 1350 
cm^^, is absent in most of our devices, as it is in Dl. This 
attests to the low-defect nature of our devices. 

We take spatially resolved Raman spectra of our devices at 
different bias voltages. The 532 nm laser with a diffraction- 
limited spot-size of 360 nm, affords 10-20 data points along 
the length of Dl. Figure 2(a) shows the Raman shifts of the 
G+ and G_ bands, along the length of the nanotube, for Vsd=0 
V and 1.2 V. At Vsd=^-'2. V, the bands downshift significantly 
and develop a spatial profile. The G+ profile fits to a fourth 
order polynomial, while a parabola suffices for the G_ profile. 
The Raman shifts observed at high bias are subtracted from 
the reference (Vsd=0 V) at every spatial point. This nullifies 
any contributions from the local nanotube environment to the 
downshift profile. 

Previously, we showed that the downshift of the G-|- and G_ 
bands could be interpreted in terms of increased phonon pop- 
ulations im and hence effective temperatures. This interpre- 
tation was in agreement with the Stokes/anti-Stokes intensity 
ratio observed in that work. In this work, we calibrated the G 
band downshifts with temperature in a temperature-controlled 
stage. This calibration data is shown in the inset of Fig. 2(b) 



for device Dl. The downshift is linear in temperature, with 
almost equal slopes for the two bands. The G band tempera- 
ture is obtained at each point along the nanotube by dividing 
the voltage-induced change in the Raman shift by the slope of 
the calibration line. The resulting temperature profiles for the 
two bands are shown in Fig. 2(b). 

This result constitutes the first observation of a spatial tem- 
perature profile of a single-walled nanotube under Joule heat- 
ing. Figure 2(b) also shows that the temperatures at the ends 
of the nanotube are higher than room temperature and highly 
asymmetric, indicative of asymmetric thermal contact resis- 
tances. Of ten 5 /im devices studied, 3 devices exhibited G+ 
and G_ bands that could be resolved separately with bias volt- 
age. All three devices showed the same qualitative behavior 
described above. Note that the G band Stokes/anti-Stokes in- 
tensity ratio was not resolvable with bias for any of these three 
devices. 

In an earlier work [llj, we reported that only one of G+ 
or G_ bands downshifted with bias. This observation was 
most prominent in shorter devices (< 2/im). Figure 3(inset) 
shows the Raman spectrum for such a 2 fim device, D2. On 
application of 14^=1 -2 V to this device, the broad G_ band 
(LO) downshifted by 22 cm^ ^ , while the other bands (G+ and 
another mode at ^^1545 cm"^) remained unchanged within 1 
cm^^, similar to results reported in [11|. Figure 3 shows the 
(converted) effective temperature of the G^ and G_ bands at 
Fsci=1.2 V, along the length of D2. Thus, while the G_ band 
shows a non-uniform temperature profile, the G+ and other 
band at 1545 cm^^ (not shown, to maintain claiity) show no 
modulation in space (within the error of the measurement). 

We now interpret these results. In thermal equilibrium, both 
G+ and G_ downshift equally, as observed in the temperature- 
controlled stage data of Fig. 2(b)-inset. Even in the case of 
thermal non-equilibrium, each phonon mode is likely to be 
at least at hot as the lattice. The lower temperature profiles 
shown in Figs. 2(b) and 3 are, thus, upper bounds on the 
lattice temperature profile. In Fig. 3, this lower profile co- 
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FIG. 3: (Color online) Spatial temperature profiles of G+ and G- 
bands. Inset: Raman spectrum at the center of D2. 
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TABLE I: Bias response of the G band (which of G-|- and G_ respond 
to bias voltage) versus approximate length L and G+/G_ intensity 
ratio for several devices. 



No. 


L (pm) 


Bias response 


Intensity ratio 


1 


0.5 


G+ 


0.34 


2 


0.5 


G+ 


0.07 


3 


1 


G+ 


0.45 


4 


2 


G- 


0.11 


5 


2 


G+ 


0.42 


6(D2) 


2 


G- 


1.53 


7 


2 


Both 


0.04 


8 


2 


Both 


0.09 


9 


5 


Both 


0.04 


10 (Dl) 


5 


Both 


0.12 


11 (D3) 


5 


Both 


3.99 



incides with room temperature; therefore, the observation that 
the G+ mode of D2 does not modulate with bias or space indi- 
cates that the lattice remains at room temperature. In contrast, 
the G_ mode exists at an elevated non-equilibrium effective 
temperature, as a direct result of generation by high-energy 
electrons. In long nanotube devices, these hot phonon modes 
transfer energy to the pool of thermalized phonons through 
inter-phonon scattering processes, heating up the rest of the 
lattice. 

We note that while the G_ bands of Dl and D2 are 
broad, indicative of strong electron-phonon coupling present 
in metallic nanotubes fTFl, their intensities differ strongly. 
This is likely the effect of differing chiralities for Dl and D2 
lfT9l . To study whether chirality plays an important role in our 
observations, we tabulate in Table |l] the bias response of the 
G band with length and G+/G_ intensity ratio for several de- 
vices. Table |l] confirms that length, rather than chirality, is the 
primary difference between equilibrium and non-equilibrium 
behavior. 

For a quantitative explanation, we first model the hot 
phonon profile for the simpler case of no lattice heating (Fig. 
3). In this case, thermal transport is not governed by Fourier's 
law (diffusive heat transport), since the thermalized phonons 
are not heated. Instead, one needs to understand the decay 
process of the hot phonons. Consider a nanotube of length 
L under uniform heating by an electric current /. Let g{I) be 
the hot phonon generation rate per unit length, d be the phonon 
decay length and t be the decay time. Hence, a phonon gener- 
ated at point X has a probability of reaching point xq given by 
exp{—\x — xol/d), which is a general characteristic of decay 
processes. Considering phonons arriving from both directions 
and a decay rate of ti/t, where n is the phonon population at 
any point, the continuity equation gives: 
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Plotting the solution for L=2 /im, one obtains a good match to 
the temperature profile of G_ in Fig. 3 using a decay length 
d ^ 1/im. 

We also need to understand the mechanism of power dissi- 
pation. In device D2, at 14(^=1.2 V and peak / ~5 /lA, the 



power input to the nanotube is 6 /iW, carried away mainly by 
thermal conduction along the nanotube. The electronic contri- 
bution to thermal conduction is known to be a small fraction 
(l/5th or less) of the equilibrium lattice contribution |20|. An 
upper bound to the lattice thermal conductance, Gth, of nan- 
otubes in quasi-equilibrium is given by the result for ballistic 
phonons flT\. At room temperature, Gth ^ Sr^ W/K |22|, 
where rt is the nanotube radius in meters. A weak RBM was 
observed in the Raman spectra of this nanotube at 240 cm^^, 
which corresponds to a nanotube diameter of 0.96 nm by the 
relation lorbm = 2QA/dt + 27 |I23|. This yields a thermal 
conductance of 4 nW/K for the nanotube assuming near equi- 
librium conditions. From the relation Q= Gth^T, the lattice 
cannot dissipate 6 /iW of power without heating to extremely 
high temperatures. However, such lattice heating is not sup- 
ported by our observation, as discussed earlier. 

We explain this apparent mismatch between heat generation 
and dissipation by considering the role of hot phonons. The 
estimate for Gth assumes that all phonon modes are in thermal 
equilibrium. Using the Landauer model for phonon transport 
to calculate the thermal power conducted EOl . 
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we can estimate an upper bound for the heat transport by hot 
modes in the ballistic phonon limit. In equation (2), ?/ is the 
Bose-Einstein distribution function, ( is the transmission set 
to 1 (see e.g. ref. |20] l, Vm is the group velocity of the mth 
phonon mode and huj^n is the phonon energy. Since the ef- 
fective temperatures under consideration (That ~900 K) are 
much smaller than hujm (~2000 K), the distribution can be 
taken to be constant, rj ~0. 1. Integrating over the range of en- 
ergies of these phonons (1300-1600 cm^^), equation (2) re- 
duces to Qph— 0-4:Nr] /iW, where N is the total number of 
hot phonon modes. Thus, it is possible to dissipate the gener- 
ated heat if ~ 150 phonon modes participate in heat transport. 
This is smaller than the total number of phonon branches of 
D2 and hence reasonable. Figure 4(inset) shows a schematic 
of this heat removal mechanism. 
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FIG. 4: (Color online) Calculated spatial variation of thermal con- 
ductivity for device D3. Inset: Proposed mechanism for thermal 
transport in short (top) and long (bottom) devices. 
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With this understanding of hot-phonon decay and thermal 
transport for the strongly non-equilibrium case, we now con- 
sider longer devices, which show conditions closer to equi- 
librium. As Fig. 4(inset) shows, the bulk of the lattice be- 
comes heated because the devices are long enough for the hot 
phonons to thermalize. In this case, the lattice could carry 
away a significant proportion of the heat generated. Since 
the lattice is almost as hot as the non-equilibrium phonons, 
we can estimate the thermal conductivity in the limit where 
all the heat is dissipated through the lattice. We use an iter- 
ative Fourier's law approach along with the Landauer model 
for electron transport, as developed in ref. [6l Following their 
formulation: 



ax ax 



dx 



(3) 



with 



Xeff{Tac{x),Top{x)) = {Kc + Kp,erns{x) + Kp,abA^)) ^ 

(4) 

Here, A is the cross-sectional area of the nanotube, k{x) 
is the thermal conductivity and R{x) is the resistance, cor- 
rected for contact resistance ^25 kfi and dependent on the 
effective electron mean-free-path Xef / (x) . Xef f{x) is a func- 
tion of the lattice temperature Tac, hot phonon temperature 
Top (which are known from our measurement) and scatter- 
ing lengths Xac.RT, Xop.min fS). The spatial variation of the 
parameters necessitates the use of a finite element iterative 
calculation. Using suitable values for the scattering lengths 
(as obtained from experiment ||2], l3l), Xac,R.T ~1.6 iim and 
Xop,mm ~180 nm, we can obtain a spatial variation of the 
thermal conductivity. Note that the thermal conductivity as 
obtained from Fourier's law (eqn. 3) is extremely sensitive to 
the diameter of the nanotube. While device Dl did not show 
an RBM in its Raman spectrum, similar data was obtained 
from device D3 (L=4.6 iim, u!fiBM=^'2-4- cm^^, dt=2.l nm). 



The above analysis yields a spatial variation of thermal con- 
ductivity for D3 as shown in Fig. 4, roughly following the 
spatial temperature profile. 

We note that the magnitudes of k are similar to other mea- 
surements on single-walled nanotubes |24|, thus validating 
the assumption in longer devices that the lattice is the dom- 
inant heat carrier. The hot-phonon thermalization length is 
thus larger than 1 /im and ^2.3 /im, considering L/2. Note 
that K shown in Fig. 4 approaches 6000 W/m-K at the max- 
imum (where T ~900 K) and is among the highest reported 
for any material |25]. 

Our data also allows a direct estimate of thermal contact 
resistance, Rth- Previously, laser heating of nanotubes was 
used to determine of the ratio of thermal contact resistances 
the left (L) and right (R) lead, Rth,Ll Rth.R |26|. In our case, 
since both temperature and heat flow are known, we can di- 
rectly compute both Rth,L and Rth.R- obtain 
Rth,L ~ 8 X 10^ KTW and RthM^ ~ 10^ K/W respectively. 
Similar values are obtained for other devices. The asymmetry 
is device specific and likely dependent on the nanotube-metal 
contact interface. Thermal contact resistance, thus, accounts 
for a significant temperature drop at the ends of the nanotube. 

In summary, the spatial temperature profile of electrically- 
heated single-walled nanotubes is obtained for the first time. 
This measurement provides insights into the mechanism of 
thermal transport and gives a spatial measure of the thermal 
conductivity and thermal contact resistances of carbon nan- 
otubes. 
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